Properties of atomic intercalated carbon K 4 crystals 
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The stability of atomic intercalated carbon K4 crystals, XC2 (X=H, Li, Be, B, C, N, O, F, Na, 
Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, Br, Rb or Sr) is evaluated by geometry optimization and 
frozen phonon analysis based on first principles calculations. Although C K4 is unstable, NaC2 and 
MgC2 are found to be stable. It is shown that NaC2 and MgC2 are metallic and semi conducting, 
respectively. 
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Recently, mathematical analysis elucidated that 120° 
bonding forms a 3-D 10-member ring periodic structure 
{K± structure) and can be regarded as the twin of the di- 
amond structure 1 -. From the view point of carbon chem- 
istry, 120° bonding is of the sp 2 bond type, which form 
a 2-D structure with a 6-member ring, and the K4 car- 
bon structure has not been experimentally discovered. 
Because of the wide applications of carbon materials 
(graphite, amorphous, diamond, CNT and fullcrene 2 -), 
if such a new member exists in nature or can be synthe- 
sized, then this new material should have a very signifi- 
cant impact. 

The stability and properties of this carbon K4 
structure have been studied using first principles 
calculations 3 -^. It is interesting to note that carbon K4 
is metallic* 3 -^, although the phonon calculation results in- 
dicate the structure is unstable^. Compared with the 
graphite structure of sp 2 bonding, the binding energy per 
bond is slightly lower. This is possibly one of the reasons 
for the instability of the carbon K4 crystals against the 
thermal vibration. 

It is well understood that atomic intercalation to the 
carbonaceous materials changes the structural stability 
and sometimes provides totally different properties. We 
propose that atomic intercalation will affect the elec- 
tronic distribution in the K4 crystals and possibly sta- 
bilize the structure. Over 40 years ago, disilicides such 
as SrSi2 were found to have K4 type flame of Si£. This 
result encouraged us to investigate impurity intercalated 
carbon systems although such crystal structures are 
not found in the carbide system^. 

This study examines the structural stability and sev- 
eral properties of the atomic intercalated carbon K4 crys- 
tals, XC 2 (X=H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, 
P, S, CI, K, Ca, Ga, Ge, As, Se, Br, Rb or Sr) by first 
principles calculations. 



First principles calculations based on density func- 
tional theory^ are performed for the K4 type crystals 
XC2, C, Si and SrSi2 using the Vienna Ab-initio Simula- 
tion Package (VASP) 1 ^. 

Local density approximation (LDA) n i 12 is used for the 
exchange-correlation energy functional. In the present 
calculations, all these crystals are considered as spin-un- 
polarizcd systems. 

To reduce the computational costs, the projector- 
augmented wave methocU^ is used to approximate elec- 
trons in each atom in the crystal. For H, Li, Be, B, C, 
N, O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, 
Se, Br, Rb and Sr, valence electrons of Is 1 , 2s 1 , 2s 2 , 
2s 2 2p\ 2s 2 2p 2 , 2s 2 2p 3 , 2s 2 2p 4 , 2s 2 2p 5 , 3s 1 , 3s 2 , 3s 2 3p\ 
3s 2 3p 2 , 3s 2 3p 3 , 3s 2 3p 4 , 3s 2 3p 5 , 3s 2 3p 6 4s 1 , 3s 2 3p 6 4s 2 , 
4s 2 4p 1 , 4s 2 4p 2 , 4s 2 4p 3 , 4s 2 4p 4 , 4s 2 4p 5 , 4s 2 4p 6 5s 1 and 
4s 2 4p 6 5s 2 are considered, respectively. 

To evaluate the stability of the XC2 K4 crystal struc- 
ture, the following procedure is adopted. 

(1) The original conventional unit cell for the carbon 
K4 crystal is shown in Figure 1(a). Four X atoms (X 
= H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, P, S, CI, 
K, Ca, Ga, Ge, As, Se, Br, Rb, or Sr) are allocated to 
reduced coordinates (|, |, ±), (|, f , §), (§, |, f), (§, 
|, |) in the original unit cell. The primitive unit cell of 
XC 2 with P4332 (O 6 ) symmetry is shown in Fig. 1(b). 
The symmetric plane of the atomic intercalated system 
is shown in Fig. 1(c). As shown in Fig. 1(d), three 
adjacent atoms (red) are aligned on a straight line with 
the same distance (C-X, and X-C). 

(2) Binding energy versus volume curve is evalu- 
ated and fitted using Murnaghan's equation of stated. 
For these calculations, Brillouin zone integration is 
performed for 8x8x8 fe-point meshes generated by 
the Monkhorst-Pack scheme. The residual minimiza- 
tion/direct inversion in the iterative subspace method is 
used to accelerate the convergence of self-consistent total 
energy calculations. The convergence criterion is set to 
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be within 1 xl0~ 8 cV/formula unit cell. The cut-off en- 
ergy for the plane- wave expansion of valence electrons for 
the primitive unit cell is determined, so that the number 
of plane waves is constant over a full range of the lattice 
constant. Around the minimum of binding energy ver- 
sus volume curve, the cut-off energies were set to be 500 
and 400 eV for C and XC2, and Si and XSi2 crystals, 
respectively. 

(3) The unit cells obtained from step (2) are optimized 
with freedom of the atomic configuration under the sym- 
metry constraint. 

(4) The unit cells obtained from step (3) are optimized 
without any constraint on the crystal structure. The con- 
vergence criterion is set to be within 1 xl0~ 7 cV/A unit 
cell. 

(5) The frozen phonon calculations were performed us- 
ing the FROPHO code^ which is based on Parlinski-Li- 
Kawazoe method^. To obtain the force constants for 
the phonon calculations, the atomic displacements are 
set to be 0.01 A . The Born - von Karman boundary 
condition is applied for each primitive unit cell obtained 
with multiplication after steps (3) and (4) for the phonon 
calculation. In these calculations, the primitive cells are 
multiplied by 2 for each direction of the unit vector. 

III. RESULTS AND DISCUSSIONS 

The structure with K4 crystal flame considered for C, 
XC2, Si and SrSi2 systems in the present study maintain 
K4 type crystal flame composed of C or Si with distortion 
through the whole process of the geometry optimization. 
Furthermore, the crystal symmetry is maintained within 
an accuracy of less than O.OOlA . The lattice constant 
(a), volume at lattice constant (Vo), cohesive energy per 
atom (E co h) and bulk modulus at Vq (Bq) are evaluated 
for the fully optimized structure. In this study, E co h for 
XC2: E co h t xc 2 i s defined as 

p -E X C. 2 +Nxato m XExato m +N Catom/Ki XEcatarn 



Here, Ea and Na are the total energy and the number of 
atoms for A, respectively. To obtain these values, Mur- 
naghan's equation of state is used to fit the binding en- 
ergy versus volume curve. Parameters in the equation are 
determined using the least squares method in the range 
of about O.8V0 < V < 1.2Vq, where the root mean square 
is set to be less than 5.0 meV/atom. 

The obtained values of a, Vo, Ecoh, and Bq for the 
crystals are shown in Figure 2(a). In XC2, the values of 
a and Vo generally increase with the atomic number of X. 
Generally, the values of a for the XC2 crystals are larger 
than that for carbon K4. That is, generally the carbon 
K4 flames are expanded by intercalating X to the carbon 
systems. 

However, the values of Vo for X=H, Li, S, are 
smaller than that of carbon K4 while those for X=C1, 
Sr arc larger. This is because of the sufficiently wide 



vacant spaces in carbon K4 that is compensated by X 
atoms and the increase of the occupied spaces by X atoms 
for the increase of the atomic number of X. 

On the other hand, the values of E co h decrease with 
the increasing period of X. There is a strong negative 
correlation with a and Vq as shown in the correlation co- 
efficients: r a _£; coh =-0.792 and ry _£ coh =-0.796. For each 
period of X in the XC2, E co h generally shows largest val- 
ues for intermediate elements X: IIIB or IVB. This en- 
hanced energy gain can be explained by the completion of 
the electronic shell in C atoms composing K4 type flame, 
the bonds formation with X in the XC2 crystals and the 
bonds cxpantion which weaken the bonds. The E co h val- 
ues of the XC2 are not larger than that of the carbon 
K4. This is apparently resulted from the relative weak- 
ness of the X-C bonds in the XC2 crystals compared with 
the C-C bonds in the carbon K4 crystal. Subsequently, 
the absence of correlation between E co h and dynamical 
stability is shown. 

The values of Bq also decrease with the increasing pe- 
riod of X, similar to the case of E co h - The Bq,s of XC2 
are generally smaller than that of C. There is a strong 
correlation between Bq and E co h as indicated by the cor- 
relation coefficient: rE aoh -B =0.848. 

From the values of a and Vq in Si and SrSi2, the relation 
shown in the C and XC2 systems also seems to hold for 
the Si and XSi2 systems. However, the values of E co h and 
Bq of SrSi2 are larger than those of Si. This relation is 
apparently different from the cases of C and XC2 systems. 
This can be attributed to the X-Si bond strength being 
much higher than that of the Si-Si bonds in XSi2- 

To understand the correlation between the stability 
and the structure in detail, the nearest-neighbour dis- 
tances between different atoms specified in Fig. 1(e) d 
C-C, d X-CU, d X-Cb, d X-Cc, d X-C (the average of all 
d X-C) and d X-X as well as the angle Z X-Ca-Cb, and 
the minimal dihedral angle for the nearest-neighbour C 
atoms Z C-C-C-C are investigated. 

As shown in Fig. 2(b), in general d C-C, d X-C and 
d X-X increase with the increasing atomic number in 
X, similar to the cases of a and Vq. The d C-C and 
d X-C show very strong correlation with the a, indi- 
cated by the correlation coefficients: r a -d C -a =0-930 and 
r -; =0.987. These distances have also correlation 

a—dx-c 

with E coh : r Ecoh -d c _ a =-0.630 and '" Bcoh _ 3jt _ c =-0.741, 
although the values are smaller. 

It is apparent that the correlations shown in 
d X-C with a and E co h are resulted from the 
very strong correlation of d X-C a and d X-Cb 
with a and E co h as indicated by the correla- 
tion coefficients: r a -d x _c =0-824, r a -d x -c B =0.973, 
r a -d x -c c =0.180, r Eaoh -dx-c A =-0.759, r Ecoh -d x -c B =" 
0.796, and r E „ k -d X - ao =-0.051. 

Here, it should be noted that the each inequivalent 
atomic distance between X and C: d X-C.4, d X-Cs, and d 
X-Cc show interesting X dependence. The XC2 systems 
with X of 3rd, 4(5)th periods show distinct differences 
between X=IA and IIA, and other Xs. In the compound 
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with X=IA or IIA, the d X-CU and d X-Cb are relatively 
larger while d X-Cc is smaller compared with those of 
the other compounds with X of the neighbour atomic 
number. The angle Z X-C^-Cs show distinct smaller 
values for the compounds with X=IA and IIA and this 
is attributed to the X dependence of the d X-CU and d 
X-C B - The dihedral angle Z C-C-C-C shown in Tab. II 
also show smaller values for the compounds with X=IA 
or IIA although the differences with other compounds 
are very small. Apparently, those differences are due to 
the number of the outer most valence electrons in X. 
However, an exceptional trend is shown in the systems 
with X of 2nd period for those values and the trend can 
be attributed to the elemental closeness between the X 
and C. 

The relative position of X in the crystal is maintained 
resulting in r a _d JC _ x =1.000. Therefore, d X-X can be 
considered as a standard distance in the XC2 systems. 
To emphasize the similarity of the structures, various 
distance ratios of the nearest-neighbour distances to d 
X-X are shown in Fig. 2(c). Although significant larger 
(smaller) d C-C / d X-X, d X-C / d X-X, and d C-C c / 
d X-X (d C-Ca I d X-X) are shown in the crystals with 
X=C, N, and O, these ratios are almost constant in the 
same group of X. 

In Figure 3, the phonon density of statcs(DOS) with- 
out step (4) for C and XC 2 (X=H, Li, Be, B, C, N, O, 
F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, Br, 
Rb or Sr), Si and SrSi2 with K4 type flame of C or Si are 
depicted. From the above discussions and Fig. 2(a), the 
intermediate value of a realized for the compounds with 
X in 3rd (4th) (4.5-4.8 A ) and the accompanied inter- 
atomic distances in these crystals seem to make the band 
narrower. Namely, the inter atomic distances roughly de- 
termine the vibrational frequency range in these crystals. 
The imaginary frequency modes rarely appear in NaC2, 
MgC2 and SrSi2 although they appear in relatively wide 
imaginary frequency range in C and Si. 

For the structures of C, NaC 2 , MgC 2 , Si and SrSi 2 
obtained after step (4), the phonon DOS and dispersion 
relationship are also evaluated. As shown in Fig. 4, the 
phonon DOS for the fully optimized structures of the 
other XC2 crystals show almost the same shapes as those 
of the optimized structures from step (3). Therefore, it 
is expected that the discussion in the previous paragraph 
is significant and skipping step (4) is expected to be ef- 
fective for reducing the computational cost in this study. 

As shown in Figure 4, in carbon K4, although the imag- 
inary frequency mode is not appear at the T point, it 
appear in wide range of other k points as indicated by 
Yao et alJ±. On the other hand, silicon K4 has imaginary 
modes even at the T point. Namely, the carbon and sil- 
icon K4 crystal structures are on a high rank of saddle 
points in those potential energy surfaces and the lifes of 
them are expected to be very short in nature. 

On the other hand, in our calculations although NaC2, 
MgC2 and SrSi2 have imaginary frequencies in the acous- 
tic modes slightly around the T point, there are no imag- 



inary modes in the other k points. These suggest that 
those structures are stable although they break the struc- 
tures against the phonon vibration with long wave length 
limit. The existence of the SrSi2 K4 in nature^ sug- 
gests that the results in this study based on the density 
functional theory can predict the thermal stability of the 
structure. From this fact and the calculation results, the 
stable existence of K4 type crystal structures in NaC2 
and MgC2 is expected in nature although they do not 
realize maximum E co h in the considered XC2 crystals. 

From the maximum frequencies obtained for diamond 
(1300 cm" 1 for G±£, 500 cur 1 for Si±2) and graphite 
(1600 cm" 1 for C)M, the orders of both C-C and Si-Si 
chemical bonds can be considered less than 1, as dis- 
cussed by Yao et a&. Therefore, based on the rough 
expectation from the mono atomic systems, the orders 
are expected to be less than 1 between the atoms com- 
posing K4 type flame of NaC2, MgC2 and SrSi2 crystals. 
The characters of chemical bonds arc analysed from the 
valence charge distribution as described below. 

To understand the stability of the K4 systems, the va- 
lence charge distributions are analysed. Fig. 5 shows var- 
ious types of distribution in the conventional or the prim- 
itive unit cells for the fully optimized C, NaC2, MgC2, 
Si and SrSi2 crystals with K4 type flame. Fig. 5(a) de- 
tails planes selected to show the distribution of quantities 
noted in Figs. 5(c), (e), and (g). The populations around 
the Sr atoms are neglected for easier comparison in Figs. 
5(c) and (e). 

Figures 5(b) and (c) show the isosurface and contour 
of charge density. Apparently, the valence charge is lo- 
cally distributed along the lines between adjacent C or 
Si atoms of the K4 type flame for the considered sys- 
tems. It should be noted that a distinctive difference ap- 
pears between the charge distribution of the C(XC2) and 
Si(SrSi2) crystals in the intermediate region between the 
nearest- neighbour atoms in those flames. The charge is 
more localized to the C than Si atoms. These features are 
also indicated in the diamond crystal structures as shown 
in the previous studies^. It can be said that the atomic 
intercalation change the distribution slightly. However, 
as discussed for the phonon frequencies, this quantity is 
sufficient to change the properties of the crystals. 

Figures 5(d) and (e) shows the differences between the 
valence charge density of the crystals (ex. pjvaC 2 ( r )) an d 
their separated components: the intercalated atoms (ex. 
PAfa(r)) and the flame (ex. pc 2 (r)). As shown in Fig. 
5(d), the significant excess charge is shown in the region 
between the adjacent Na-C, Mg-C and Sr-Si for NaC2, 
MgC2 and SrSi2 , respectively; therefore suggesting bond 
creation in those regions. Furthermore, as shown in Fig. 
5(e), in NaC2, MgC2 and SrSi2, the charge accumula- 
tions increase in the region of n type orbitals while those 
decrease along the a type bonds; suggesting that the 7r 
type bonds are strengthened and the a type bonds are 
weakened. On the other hand, the maximum frequency 
of the phonon modes shown in Fig. 4 is less than the 
mono atomic systems. Therefore, the net bond orders in 



4 



the K4 type flame are expected to be decreased by the 
atomic intercalation. 

Here, the nature of the bonds is noted from the infor- 
mation of the ELF. Figs. 5(f) and (g) show the isosur- 
face (0.8 [-]) and the contour. As shown here, C, NaC2, 
MgC2, Si and SrSi2 has similar ELF distribution between 
the adjacent C-C and Si-Si. However, the feature of ir 
type bonds formation is strengthened with significance. 
The relatively smaller value of Z C-C-C-C (Z Si-Si-Si-Si) 
seems to have correlation with this feature. The strength- 
ened 7r type bonds seem to stabilize the structures against 
the themal vibration. From the comparison with the case 
of graphite, the absense of the tt type bonds indicated by 
Yao et aL— , is not confirmed in our calculations. 

The stabilization of atomic intercalated carbon crystals 
is realized from delicate balance of bonding. As shown in 
the valence charge density, the formation of strongly po- 
larized a type bonds between Na or Mg and Ca and the 
strength enhanced tt bonds seems to stabilize the struc- 
ture of NaC2 and MgC2- The discriminative stability 
seems to be realized by the properous atomic radii of X 
and the completion or near completion of the electronic 
shell of the each C atom in the K4 type flame by X=IA 
or IIA. Similar mechanism of the stability is expected to 
exist in the Si and XSi2 systems. 

Figure 6 shows X-ray diffraction (XRD) patterns of the 
fully optimized crystal structure of C, NaC2, MgC2, Si 
and SrSi2- Monochromatic radiation with wave length 
1.541 A is assumed in these calculations. As shown in 
this figure, the compounds have specific XRD peaks at 
larger 26 compared with the pure K4 crystals. 

Figure 7 shows the electron band structures, DOS, and 
local DOS of valence electrons for the atoms composing 
K4 type flame and the intercalated atoms in the fully op- 
timized C, NaC2 and MgC2 and the Si and SrSi2 crystals, 
respectively. From the bottom energy levels to the upper 
levels, the angular momentum for the electronic states 
successively changes in the order of s, p and d character, 
although every component appears in most states. From 
the bottom energy levels to those around the Fermi level, 
the band structures are similar in these crystals. NaC2 
and MgC2 K4 show a metallic and semi conducting fea- 
ture, respectively, while carbon K4 is metallic. As shown 
in Fig. 6, both Si and SrSi2 have metallic features. It is 
noteworthy that the previous first principles calculations 
for SrSiz^ also show qualitatively the same result with 
the present calculation. 



IV. CONCLUSIONS AND SUMMARY 

The stability of the atomic intercalated carbon K4 
crystals, XC 2 (X=H, Li, Be, B, C, N, O, F, Na, Mg, 
Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, Br, Rb, and 
Sr), are evaluated from the geometry optimization and 
the frozen phonon analysis based on the first principles 
calculations. NaC2 and MgC2 are found to be stable al- 
though C K4 is not stable. The formation of strongly 



polarized a type bonds between Na or Mg and Ga and 
the strength enhanced tt bonds seems to stabilize the 
structure of NaC2 and MgC2 . The discriminative stabil- 
ity seems to be realized by the properous atomic radii of 
X and the completion or near completion of the electronic 
shell of the each C atom in the K4 type flame by X=IA 
or IIA. The NaC2 and MgC2 are found to be metallic 
and semi conducting, respectively. 

V. APPENDIX 

The values presented in Figure 2 are listed in Table I 
and II. 
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FIG. 1: (a) The conventional unit cell of the carbon K4 crystal with I4i32(0 ) symmetry, (b) Initial structure of the primitive 
unit cell of XC2 crystal with K4 flame of C with PAsZ2{0 6 ) symmetry, (c) 2x2x2 super cell of the (b). The symmetric plane 
of the impurity intercalated system, (d) Three adjacent atoms (C-X-C) that aligne in a straight line with equivalent distance 
(C-X and X-C) are encircled with red thick line in (c). (e) An extracted structure. The C atoms can be classified into Ca, Cs 
and Co- 
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FIG. 2: (a) Determined lattice constant(a) , volume at a(Vo), cohesive energy (E co h), and bulk modulus at Vo(Bq), and (b) 
the nearest-neighbour distances and angles for the fully optimized C, XC2, Si and SrSi2 crystal structures with K4 type flame 
within LDA. The classification of the C(Si) atoms (Ca, Cb, and Cc)) is based on the irreducibility of the crystal symmetry as 
shown in Fig. 1. (c) The distance ratios. 
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FIG. 3: Phonon density of states(DOS) of XC 2 (X=H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, P, S, CI, K, Ca, Ga, Ge, As, Se, 
Br, Rb and Sr), C, Si and SrSi2 crystal with K4 type flame of C, or Si. The structures are obtained from geometry optimization 
under symmetry constraint to the initial structures. 
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FIG. 4: (Left) Phonon density of states (DOS) of partially and fully optimized structures of C, NaC2, MgC2, Si and SrSi2 with 
K4 type flame of C or Si. (Right) The phonon dispersion relationship for the fully optimized structures. 
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FIG. 5: (a) The conventional or the primitive unit cells for the fully optimized C, NaC2, MgC2, Si and SrSi2 crystals with K4 
type of flame of C or Si. The planes are selected to show the distribution of the valence charge density, difference between the 
valence charge density of the crystals (ex. pjvaC 2 ( r )) an d their separated components: the intercalated atoms (ex. pjvo(r)) 
and the flame (ex. pc 2 ( r ))> an d the electronic localization function (ELF). The isosurfaces of the valence charge density (1.5, 
1.0, 1.0, 0.35 and 0.35 e/A 3 for C, NaC2, MgC2, Si and SrSi2, respectively) (b) and the contours (c). The isosurfaces of the 
differences of the valence charge density (0.05, 0.05 and 0.025 +(-)e/A' ! coloured with red(blue) for NaC2, MgC2 and SrSi2, 
respectively) (d) and the contours (e). The isosurfaces of the ELF (0.8 [-]) (f) and the contours (g). For (c) and (e), the high 
population around the Sr atoms is neglected for easier comparison. XCrySDen^i was used for the visualization. 
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FIG. 6: XRD patterns 



for the fully optimized structures of C, 



NaC 2 , MgC 2 , Si and SrSi 2 with K 4 type flame of C or Si. 
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FIG. 7: Electronic density of states (DOS) and the dispersion relationship for the fully optimized conventional or primitive 
unit cells of C, NaC2, MgC2, Si and SrSi2 crystals with K4 type flame of C or Si. 
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TABLE I: Determined lattice constant(a), volume at a(Vo), cohesive energy (E co h), and bulk modulus at Vb(-Bo) for the fully 
optimized C, XC2, Si and SrSi2 crystal structures with K4 type flame within LDA. 



Component 


a [A] 


V [A 3 /atom] 


E co h [eV/atom] 


Bo [10 11 &] 


oi 


4.082 


8.502 


7.702 


2.67 


HC 2 


4.025 


5.433 


6.335 


3.00 


T,i(~!o 


4.250 


6.396 


6.231 


2.41 


j_)LV^2 


4.412 




fi 4D1 


9 to 


RP. 
iJv^2 


4 ^98 


fi 7^8 


7 9^0 


2.94 


ecu 


4.326 


6.747 


7.295 


3.16 




4.348 


6.850 


7.003 


3.33 




4.218 


6.252 


7.221 


O.U 1 


FC 2 


4.296 


6.607 


6.273 


2.98 


NaC2 


4.526 


7.728 


5.476 


2.06 


ivigv^2 


4.506 


7.623 


5.856 


2.29 


rtlv^2 


4 ^QS 


8 1 09 


6.214 


2.29 


01\^2 


4.622 


8 997 


fi TO1 


9 TO 


PP. 
A v_>2 


4.619 


8.214 


5.805 


2.42 




4.644 


8.345 


5.261 


2.38 


CIC2 


4.699 


8.646 


5.152 


2.18 


av»^<2 


4 Q9fi 




4 ROfi 


1 TO 


CaC 2 


4.803 


9.232 


6.099 


1.82 


GaC 2 


4.679 


8.536 


5.831 


2.41 


GeC 2 


4.731 


8.823 


5.875 


2.19 


a or 1 
ASU2 


4. M O 


y.Uoy 


0.450 


on 


SeC 2 


4.823 


9.351 


4.987 


2.10 


BrC 2 


4.908 


9.852 


4.668 


1.86 


RbC 2 


5.175 


11.550 


4.172 


0.95 


SrC 2 


5.096 


11.030 


5.545 


1.44 


Si 


6.354 


32.089 


4.220 


0.46 


SrSi 2 


6.436 


22.216 


4.634 


0.66 
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TABLE II: Determined nearest-neighbour distances and angles for the fully optimized C, XC2, Si and SrSi2 crystal structures 
with K4 type flame within LDA. The classification of the C(Si) atoms (Ca, Cb and Cc)) is based on the irreducibility of the 
crystal symmetry as shown in Fig. 1. The smallest dihedral angles for the nearest-neighbor C(Si) atoms of K4 flame are also 
shown. 



Component 


— 

d C-C [A ] 


— — — 

d X-C [A ] 


— 

d X-Ca [A ] 


— 

d X-Cb [A 


] d X-C c [A ] 


— 

d X-X [A ] 


Z X-Ca-Cb [°] 


Z C-C-C-C [°] 


C£ 


1.443 


2.209 


1.768 


2.282 


2.282 


2.500 


90.000 


70.529 


HC 2 


1.492 


2.226 


1.350 


2.180 


2.563 


2.463 


100.081 


67.946 


LiC 2 


1.525 


2.286 


2.071 


2.444 


2.199 


2.601 


84.183 


69.686 


BeC 2 


1.581 


2.409 


1.683 


2.417 


2.644 


2.701 


95.477 


69.782 


BC 2 


1.582 


2.383 


1.526 


2.354 


2.698 


2.650 


98.434 


68.732 


CC 2 


1.628 


2.405 


1.391 


2.341 


2.808 


2.649 


101.393 


67.330 


NC 2 


1.679 


2.437 


1.295 


2.347 


2.908 


2.661 


103.471 


67.202 


oc 2 


1.582 


2.343 


1.365 


2.282 


2.729 


2.582 


101.180 


65.807 


FC 2 


1.524 


2.333 


1.741 


2.372 


2.492 


2.628 


92.955 


70.313 


NaC 2 


1.614 


2.436 


2.152 


2.585 


2.381 


2.769 


85.408 


70.005 


MgC 2 


1.628 


2.422 


2.244 


2.608 


2.295 


2.757 


83.057 


69.323 


A1C 2 


1.624 


2.483 


2.006 


2.572 


2.554 


2.813 


89.592 


70.525 


SiC 2 


1.637 


2.508 


1.890 


2.555 


2.666 


2.827 


92.542 


70.369 


PC 2 


1.634 


2.503 


1.916 


2.559 


2.643 


2.825 


91.904 


70.439 


sc 2 


1.645 


2.520 


1.920 


2.573 


2.666 


2.843 


92.100 


70.420 


C1C 2 


1.663 


2.548 


1.938 


2.601 


2.698 


2.874 


92.170 


70.412 


KC 2 


1.752 


2.649 


2.322 


2.804 


2.602 


3.011 


85.787 


70.088 


CaC 2 


1.730 


2.584 


2.369 


2.773 


2.466 


2.940 


83.574 


69.498 




l.OOO 


z.ooU 


z.uzu 


O CIO 


Z.Ol ( 


z.ooo 


on 1 ni 


/ U.ozy 


GeC 2 


1.673 


2.556 


2.082 


2.653 


2.617 


2.896 


89.211 


70.513 


AsC 2 


1.685 


2.578 


2.071 


2.667 


2.659 


2.919 


89.841 


70.528 


SeC 2 


1.704 


2.608 


2.069 


2.689 


2.707 


2.951 


90.387 


70.525 


BrC 2 


1.736 


2.659 


2.067 


2.728 


2.787 


3.004 


91.257 


70.490 


RbC 2 


1.845 


2.786 


2.453 


2.954 


2.729 


3.167 


85.575 


70.042 


SrC 2 


1.834 


2.740 


2.507 


2.939 


2.619 


3.118 


83.684 


69.533 


Component 


d Si-Si [A ] 


d X-Si [A ] 


d X-SU [A ] 


dX-Sis [A 


] d X-Sic [A ] 


d X-X [A ] 


Z X-Si. 4 -Sis [°] 


Z Si-Si-Si-Si [°] 


Si 


2.247 


3.439 


2.752 


3.553 


3.553 


3.891 


90.000 


70.529 


SrSi 2 


3.463 


2.358 


3.324 


3.772 


3.200 


3.941 


81.270 


68.605 



